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Abstract 

In recent studies, the production rate of photons or lepton pairs by 
a quark gluon plasma has been found to be enhanced due to collinear 
singularities. This enhancement pattern is very dependent on rather strict 
collinearity conditions between the photon and the quark momenta. It was 
estimated by neglecting the collisional width of quasi-particles. In this 
paper, we study the modifications of this collinear enhancement when we 
take into account the possibility for the quarks to have a finite mean free 
path. Assuming a mean free path of order (g 2 T ln(l/<?)) _1 , we find that 
only low invariant mass photons are affected. The region where collision 
effects are important can be interpreted as the region where the Landau- 
Pomeranchuk-Migdal effect plays a role in thermal photon production 
by bremsstrahlung. It is found that this effect modifies the spectrum of 
very energetic photons as well. Based on these results and on a previous 
work on infrared singularities, we end this paper by a reasonable physical 
picture for photon production by a quark gluon plasma, that should be 
useful to set directions for future technical developments. 



LAPTH-790/2000, BNL-NT-00/11 



1 Introduction 



The photon production rate is thought to be a quantity of phenomenological 
interest in heavy ions collisions, possibly enabling one to detect the formation 
of a quark-gluon plasma. Part of the interest in this electromagnetic observ- 
able comes from the fact that photons are relatively weakly coupled to nuclear 
matter (a EM <C ct s ). Given the typical size of the system in such collisions 
(much smaller than the mean free path of a photon) , they do not re- interact be- 
tween their production and their observation. As a consequence, photons (real 
photons, or virtual photons decaying eventually into a lepton pair) can provide 
information on the state of the system at the time they were produced. 

In order to calculate the photon yield from a hot quark-gluon plasma, ther- 
mal field theory is the tool of choice since its Feynman rules automatically take 
into account the presence of a thermal bath with the appropriate distributions 
of partons. Thermal gauge theories are however plagued by infrared singulari- 
ties arising from the Bose-Einstein distribution functions which are singular at 
zero energy. An improvement over the bare Feynman rules is achieved by the 
resummation of one-loop thermal contributions known as hard thermal loops 
(HTL) |jj 0, 0, These thermal corrections make it possible to include in 
the propagators effects like Debye screening or Landau damping, and transform 
partons into massive quasi-particles. From a quantitative perspective, they pro- 
vide important changes to the dynamics of soft modes (of momentum of order 
gT or less). This resummation lacks however two features: it does not provide 
any Debye screening for static magnetic fields (such a screening is expected to 
arise non-perturbatively in QCD at the length scale (g 2 T)~ 1 ), and its quasi- 
particles do not undergo collisions (their collisionalQ width is also a sub-leading 
effect of order g 2 T \n{l / g)) . 

In the framework of thermal field theory, the photon/dilepton rate is ob- 
tained via the calculation of the imaginary part of the retarded photon po- 
larization tensor H, fjj. This object has been evaluated at one-loop in the 
HTL-improved perturbative expansion, both for virtual ^| and real pho- 
tons Jl^, H), [l3| Q . More recently, new processes like bremsstrahlung were 
studied in detail in this framework and have been found to be dominant sources 
of low invariant mass photons Jl^, [l6|, 0. Despite the fact that this process 
arrives only in 2-loop contributions to the photon polarization tensor, it is al- 
ways important because of a strong collinear enhancement. Indeed, it was found 
in fill that for a very small photon invariant mass, the corresponding diagram 
contains collinear singularities that, when regularized by a thermal quark mass 
of order gT, give an enhancement by a factor of order 1 /g 2 over naive estimates 
coming from power counting. 

'In a plasma, it is very important to distinguish two contributions to the total width: the 
decay width made of the zero temperature contribution to the imaginary part of the quark 
self-energy, and the collisional width which exists only in a medium. The latter is also called 
"anomalous damping rate" (or just "damping rate" for short) and is of order g 2 T\n(l/g), 
while in QED/QCD the former starts at the order g 4 . The width T that we introduce in this 
paper is the collisional width, and its inverse is the mean free path of the quark in the plasma. 
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After finding that certain 2-loop terms are contributing at leading order, one 
may wonder whether this result is specific to this 2-loop contribution only or if, 
on the contrary, this is an indication of the breakdown of perturbative expansion 
(even improved with HTLs) . In a recent paper |T^] , we studied what is the effect 
of loop corrections to this 2-loop diagram. Power counting indeed indicates a 
problem very similar to the problem raised by Linde for the calculation of 
the free energy, due to the lack of Debye screening for static magnetic modes. 
This problem is avoided for the production of massive enough photons, because 
some cancellations (occuring within any given topology, when one is summing 
over all the cuts contributing to the imaginary part) generate a kinematical 
cutoff large enough to prevent any sensitivity to the non-perturbative scale g 2 T. 
Unfortunately, this cutoff is smaller than g 2 T whenever the invariant mass of 
the produced photon is too small (typically Q 2 < g 2 Tqo for qo < T). In this low 
invariant mass region, the photon rate is therefore non perturbative: exchanged 
transverse gluons reach the scale g 2 T of the non-perturbative "magnetic mass" , 
and an infinite set of diagrams must be resummed. 

In this paper, we present a completely different approach to higher order 
corrections, that completes the picture outlined in [l8| . The idea behind the 
present study is that a width on the quark propagator will act as a regulator 
in the collinear sector, because it moves the poles of the propagator away from 
the real energy axis. Such a collisional width is necessarily a higher-loop effect, 
because the hard thermal loop framework does not take into account the colli- 
sions of quasi-particles. Having in mind the fact that 2-loop contributions are 
important because of collinear enhancement, an important question to answer is 
how much of this enhancement is lost when an additional regulator like a width 
is taken into account. This is the question we want to address in this paper, by 
calculating the same 2-loop diagrams as in jl5[ in the presence of a quark 
width. To be more definite, and keep the model as well as the calculations 
simple, we use a momentum-independent width. 

A word of caution is necessary here: the formulae found using this simple 
model should not be taken as an accurate quantitative account of what the 
effect of such a width will be on thermal photon production rates. Indeed, the 
constant width model is probably too naive to be realistic, and more importantly 
our calculation disregards the fact that a modification of the vertices should in 
principle accompany the modification of the quark propagators. Nevertheless, 
this simple approach is sufficient here for our purpose which is just to determine 
the region in which effects of a width of order g 2 Tln(l/g) are to be expected, 
since that gives another handle on how and when higher order corrections may 
affect photon production by a quark-gluon plasma. 

Wc find that the region where a width of order g 2 T\xi(\/g) is important is 
very similar to the one found in |l8| ] for the contribution of higher-loop topolo- 
gies due to an IR sensitivity to the scale of the magnetic mass. Despite their 
similarity, these two non perturbative regions have different physical interpre- 
tations. In particular, we find that the sensitivity to the collisional width of the 
quarks is a manifestation of the Landau-Pomeranchuk-Migdal effect. Indeed, it 
occurs when the formation time of the photon is larger than the mean free path 
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of the quark producing the photon. An interesting consequence of our study is 
that the LPM effect also modifies the spectrum of highly energetic photons. 

The structure of the paper is as follows. Section |^ makes more precise 
our modeling of the mean free path for the quarks. In section we start by 
computing the 1-loop contribution in presence of a width. Although not related 
to collinear singularities in any way, the purpose of this warmup exercise is two- 
fold: it illustrates the technology (how one does calculations with a width), and 
it shows in a simple way how collisions can open up the phase space. We also 
show that this 1-loop contribution is canceled by the resummation of vertex 
corrections, and is therefore not physical. 

In section ^, we present the 2-loop calculation with a quark width, and 
obtain a rather simple generalization of the formulae of flfifl . These formulae 
are discussed extensively in section |^, in which we also determine the region 
where a width of order g 2 T\n.(l/g) plays an important role. 

Section ^| is devoted to establishing the connection between the previous 
results and the LPM effect. In section |?j, we study the process obtained from 
bremsstrahlung by crossing symmetry, which turns out to be important in the 
region of large photon energy. We show that this process is also affected by the 
LPM effect. The last section contains concluding remarks. In particular, we 
combine the present work on the LPM effect with previous results on infrared 
singularities Q in order to make a syntesis and present a reasonable physical 
picture of thermal photon production. 



2 Model 

Let us first make our framework more definite. The only modification compared 
to the context extensively described in fl6|| is that the quark propagators are 
given a width, as outlined by the following substitution in the retarded and 
advanced propagators: 

Ar < a {P) = P*-M* o ±i P0 e ~" Ar - a {P) = {p Q ±iT) 2 -p 2 ~M 2 O0 ' (1) 

where Moo ~ gT is the usual asymptotic thermal mass for hard quarks |2^, [l6| , 
and where T is a constant width. The expression in Eq. (|l|) is sufficient for the 
physics we consider, which is dominated by hard quark momenta close to the 
mass shell. Whenever we need to estimate the order of magnitude of a term 
containing T in the following, we assume that it is of order g 2 T\n{\/g) (with 
g <C 1). Note that we could as well have written 

Ar - a (P) = P 2 - Ml ± 2i Po T (2) 

since the two differ by the small correction T 2 -C M<^ to the real part of the 
denominator. The retarded (resp. advanced) propagator has two complex poles 
in the po plane, located at po = ±u p — iT (resp. po = ±o-> p + iT), with u> p = 
V(P 2 + Ml). 
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3 1-loop study 



3.1 Calculation 

Let us now consider the simplest calculation conceivable in this framework, 
namely the computation of the 1-loop contribution to the photon polarization 
tensor. Our purpose is simply to illustrate possible changes brought in when 
taking into account a width for the quarks running in the loop. We do not at- 
tempt a complete calculation including HTL vertices, as pictured on Fig.g. The 
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Figure 1: 1-loop diagram. The boldface quark line denotes the resummation of 
a width on the fermion propagator. 

contribution of this diagram to the imaginary part of the retarded polarization 
tensor is trivial to obtain: 

ImIlVQ)« le 2 J jLL[n F (r )-n p ( Po )} Tr 

x[A R (P)- A A (P)] [A R (R) -A A (R)}, (3) 
where Tr denotes the Dirac's trace associated to the quark loop: 

Tr = -8P ■ R. (4) 

The case with T = is made simple by the fact that the differences A R — A A 
are d functions, thereby enabling some of the integrals to be performed trivially. 
Let us remind that for the case Q 2 < 4M^ which we assume throughout this 
paper, ImlT 1 ^(Q) = if T = due to incompatibilities between the S functions. 
To evaluate Eq. (^) when r ^ 0, we perform the integral over pg by closing the 
real axis in the complex energy plane and using the theorem of residues. 

Note first that one can disregard the poles of the statistical weights n F (po) — 
n F (ro). Indeed, these poles are the imaginary fermionic Matsubara frequencies 
and are to be plugged into differences like A R (P) — A A (P). The fact that 
r«T makes these differences very smallQ. In other words, the only important 

2 Strictly speaking, these contributions are needed to ensure that the final result is a real 
number. Indeed, when one picks poles like po = + *r from the propagators and plugs them 
in the distribution functions, the latter become complex numbers. Their (small) imaginary 
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terms are those for which denominators like P 2 — M% a are small, and comparable 
to pcT. This cannot happen if po is an imaginary number of order T. 

We are therefore left with the poles of the propagators themselves. At this 
point, we obtain the following result: 



1 f d?p - 

ImlP^Q) w -e 2 J ^3 2J K(?o + W P ) " M^p)] Tr(p = r)W p ) 

(5) 



?7 

x— — 

2cj p 



»7=±1 

1 1 



(qo + ?7W P + 2iT) 2 - uj 2 (q + r]u p - 2iT) 



We can also note for later use that the result of this approximation could have 
been obtained by starting from an expression slightly different from Eq. (^): 

Imn%(Q) « \e 2 J i-^ [n F (ro) - n F (p )] Tr 

x2^( Po )5(P 2 - Mi) [/A 2r (i?) - (R)] , (6) 

where we denote 

Ar - a ^ = (r ± 2iT) 2 -r 2 - Mg ' ^ 

In other words, if T -C T, then we can as well put twice the width on one of the 
two quark lines, and nothing in the other quark line^]. 

After Eq. (|J), the angular integration over dfl p is trivial. In the case where 
r\ = +1 (process q — > 57), we find: 

+00 

e 2 r /" 

Im 11% (Q) ps — / rfp(5 + p) [ra F (q + p) - n F (p) 



xln( (. P + p)V + r 2 ( go+ ^ 
V(wp-p)V + r 2 (©+p)V 

In particular, assuming for the sake of simplicity that q ps qq, it is trivial to 
obtain the following asymptotic behaviors for soft and hard photons: 

If q «T, Imn%(Q)^e 2 rTln(l + 4 (7 2 /r 2 ) , (9) 
If 9o>7\ ImlP^Q) — e 2 rTln(T 2 /T 2 ) . (10) 

The main point is that these contributions are proportional to the width and 
vanish when T = 0, in agreement with a direct calculation. As a side remark, 
one may note that for T ~ g 2 Tln(l/g) and soft photons (qo ~ gT), Eq. (||) is 

part is canceled by the (small) contribution coming from the poles of n F (ro)— n F (po)- For this 
approximation to be consistent, one must also neglect iF in the hard argument of statistical 
weights. 

3 This is where it is important to have a constant width. Indeed, some intermediate step 
relies on the cancellation T(P) — = 0. 
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larger by a factor 1/g than the 1-loop HTL result, while in the regime of Eq. ( flO| ) 
it is of the same order. It is the collision partners of the quarks that open up the 
phase-space (hard quarks colliding in the medium can emit a photon, a process 
forbidden for collisionless quarks) and make these contributions so large. 

Another property of this result is that there is a suppression if qo <C T, such 
that fonW* ^(Q) tends to when qq — > 0. Equations similar to Eqs. (0) and ( pl)| ) 
can be obtained when considering the contribution rj = — 1 to Eq. m). 

3.2 Cancellation with vertex corrections 

This suppression was interpreted in [5lJ as the manifestation of the LPM effect. 
However, despite the suppression at small qo, the connection with the LPM effect 
is not clear in this context. Indeed, the LPM effect is expected when the photon 
formation time is larger than the quark mean free path, a condition which never 
appears in H), nor in the above 1-loop calculation. In fact, following [^2], |23|, 
we know that the propagator of a quark close to its mass shell on which a 
width is resummed can be evaluated by an eikonal resummation of the multiple 
collisions. By the same method, one can include in this resummation all the 
vertex corrections. Indeed, one can check that the photon polarization tensor 
at this level of approximation is proportional to (in space-time coordinates) 

U^(x, y)=e 2 J [dA^e^Tr [fS(x, y\AyfS(y, x\A)} , (11) 

where S[A] is the action of the gauge fields^], and where S(x,y\A) is the prop- 
agator of a quark in the background field A. This formula includes only one 
quark loop (in addition to the quark loops that may have been resummed in the 
gluon propagators and vertices in jSf-A]), and all orders in the gluon fields. In 
the eikonal approximation, the quark propagator has a very simple dependence 
on the field A^ (which is a matrix T a A^ in QCD): 

y° 

S(x,y\A) = S (x iy )Vexp-ig J v^dt , (12) 

x° 

where is the 4- velocity of the quark, and So(x,y) is the free propagator of 
the quark. Now, if one is looking at photons of very small invariant mass, the 
two quarks are nearly parallel in the collinear limit (even if the photon is hard) , 
so that the two quarks have mostly the same v. As a consequence, the two path 
ordered exponentials cancel each other, and the product of the propagators 
under the functional integral is 

rS(x,y\A)^S( y, x\A) » 7 »S (x, y)-fS (y, x) . (13) 

4 This action does not play any role in the argument, and therefore can even include the 
HTL effective action for gluons. 
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Therefore, at the level of approximation at which the resummed quark propa- 
gator is calculated, the sum of all the gluon loop corrections is vanishing: 



all gluon corrections 




(14) 



This result generalizes a result already known for QED established by 24 |25|, 
p6| , |2?J for the subset of ladder vertex corrections, and extended to all abelian 
topologies in [p8[ , to the case where gluons are exchanged (i.e. to non abelian 
topologies). Indeed, we see that what makes gluons specific, namely the fact 
that they can couple to each other, is hidden in the action S[A], which plays a 
passive role in the argument. 

This result is also closely related to the fact that there are no jjg ■ ■ ■ g HTL 
vertices (see [^9| for an interpretation of the cancellation found in [2^, as 
a consequence of the absence of HTL vertices with n > 2 photons), because the 
sum of all the eikonal contributions can also be written as (when Q is soft): 



# \ A I A f » 

% J I f 1 I J 

TV' 1 ' [Q) = -Jw| Ww- + ^^!S»!S^v + ^vStolaS^n, + w^J»«^W +... , (15) 

t \ 



where the sum is extended to all the possible ways to close the external gluonic 
legs of the 775 • • • g vertices. If is then obvious to see that the sum reduces to 
its first term if these vertices are vanishing. 

One is therefore left with the bare 1-loop diagram, which does not contributes 
to the imaginary part of the photon polarization tensor. In other words, the 
results derived in Eqs. (|^) an d (|lC| ) as a warmup and that were claimed to be 
related to the LPM effect in |21[, are just artifacts with no physical meaning. 
In section |f] of the present paper, we show where the LPM effect appears in 
thermal field theory. 

Despite the fact that this contribution is not physical, one learns two things 
from this calculation: (1) one has to keep terms beyond the eikonal approxima- 
tion in order to circumvent this cancellation, and (2) the width may open the 
phase space to new processes. 



4 2-loop calculation 

The lesson from the previous section is that physical contributions must be 
looked for beyond the eikonal approximation. This implies inserting explicitly 
a gluon exchange with momentum L in the diagram, and not assuming that 
I? < 2P- L. 

There are in principle two 2-loop topologies contributing to the photon po- 
larization tensor. However, it was found in Jig] that only the topology correcting 
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the qq^y vertex (see Fig. is relevant in the region where the collinear enhance- 
ment takes place. Since our purpose it to study how this collinear enhancement 
is affected by the width T, we limit the present study to the terms that were 
found important in |l6fl. 



Figure 2: 2-loop vertex diagram. The dotted lines are the various cuts con- 
tributing to the imaginary part of the photon polarization tensor. The boldface 
quark line denotes the resummation of a width on the fermion propagator. The 
dot on the gluon indicates the resummation of HTLs on its propagator. 

There is though one difference with the case T = that we must take into 
account: with a zero width, processes corresponding to cuts (a) and (d) are 
kincmatically forbidden for small Q 2 . As in the previous section, turning on the 
width opens those production channels, and we cannot disregard them a priori. 
This is not a problem though from a purely technical point of view, because 
it turns out that the sum of all the cuts has a simpler expression than each 
individual pieceQ. 

If we select only terms that have the large Bose-Einstein factor n B (l ), the 
contribution of Fig. || to the imaginary part of the photon polarization tensor 
takes the following simple form Jl6| 

Imn%(Q)«ieV J AJL-[n p (r )-n p (p )} 

d A L 
(2tt 



x / 77^n B {h)p TtL [L)P^ L [L)'Tr pa 



[A R (P)A R (P + L)- A A (P)A A (P + L)] 
[A R (R)A R (R + L) - A A (R)A A (R + L)] , (16) 



where p T L (L) are the spectral functions of transverse and longitudinal gluons, 
P p<7 (L) are the corresponding projectors, and where Tr po - is the result of the 
Dirac's trace for the quark loop. The result of the previous section requires that 
we keep in this Dirac's trace only terms that do not appear in the eikonal approx- 
imation (i.e. which includes soft corrections to the hard loop momentum). In 

5 In the appendix |a[ we show how one can obtain the difference between the cuts (c) and 
(d). Having already their sum, we can therefore reconstruct the two contributions separately. 
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the collinear limit, the first non vanishing term (beyond eikonal approximation) 
is 

Tr pCT « ~&L 2 (R p R a + P p P a ) , (17) 

which turns out to be the same as the term found in |[6[ . Its contractions with 
the projectors are given by 

PP"(L)Tr pa « -8L 2 (r 2 +p 2 ){l - cos 2 9') 

PP a (L)Tr pa » +8L 2 (r 2 + p 2 )(l - cos 2 9') , (18) 

where 6*' is the angle between the 3- vectors p and I. At this point, we have used 
the fact that we are looking at collinearly enhanced terms, and consider that p 
and V are parallel (the only place where we do not do this approximation is in 
the denominators which are very sensitive to the details of the collinear sector) . 

In order now to perform the integral over po, we follow the method of the 
previous section, and use the same approximations concerning the statistical 
weights. In addition, we compute only the contribution of cuts (c) + (d), and 
multiply the result by an overall factor 2 in order to take into account the 
contribution of the cuts (a) + (b). Following the remark of the previous section, 
we can start directly from the expression: 

ImIP„(Q) » l -e 2 g 2 J ^ [n,(r„) - n,(po)] 

^n B (l )p TL (L)P^ L (L)T rpa 



(2tt 

x2ne(p )5(P 2 - Ml) 



(P + Lf-Ml 
xT>isc[Af(R)Af(R + L)} , (19) 

where we use the notation Disc/(r) = f(T) — f(—T). We first do the po inte- 
gration for free thanks to the S(P 2 — M 2 ^). In this section, we consider only the 
case of po = +ll> p (bremsstrahlung of a quark) in order to keep the calculation 
compact. The contribution of po < is identical if qo <C T, but is different if go 
is large. In the latter case, the corresponding process is a q*q annihilation (g* 
denotes a quark placed off-shell by a scattering) instead of bremsstrahlung (see 
Fig. ||), and is considered in more detail in section [?]. 

Then, it happens that the angular integral over the direction f2( of the 3- 
vector I can be done analytically in a rather simple way. We have to perform 
an integral like[] 

_fdO.il 1 
II = J 17 2L ■ P + L 2 2L ■ R + L 2 + 2P ■ Q + Q 2 + 4ir T ' (20) 



6 Strictly speaking, we have also a factor of 1 — cos 2 6' in the numerator that depends on 
f2;. This factor can be taken into account analytically in the angular integral, the price to pay 
being more cumbersome expressions. One can however make the following simplification: if 
the width is small and if we use the collinear approximation for the numerators, then this angle 
is approximately given by cos 8' si lo/l, and does not play any role in the angular integrals. 
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Figure 3: Different processes contained in the 2-loop diagram of Fig. [|. Left: 
po = +ll> p > 0. The corresponding process is the bremsstrahlung of a quark. 
Middle: po = —lu p < and lu p > qo so that vq ~ Po + Qo (and also tq + Iq) 
is negative as well. This region contains the bremsstrahlung of an antiquark. 
Right: po — —ujp < and lu p < qo so that ro = po + qo (and also tq + Iq ) is 
positive. The corresponding process is the annihilation of an antiquark with an 
off-shell quark. 



which can be rewritten as 



dCl, 



4tt 2L ■ A 2L ■ B ' 
where L = (1,1) provided we introduce the fictitious "4- vectors" 



(21) 



A= (p l + —,lp) 
L 2 

B = Mo + ^r + P-Q 



Q 2 



2ir Q T,lr) 



(22) 



The advantage of rewriting I L like this lies in the fact that the last integral is 
known in closed formQ: 



I, = 



8VA 



\n(A ■ B + VA) - ln(A ■ B - \/A) 



(23) 



where A = (A ■ B) 2 - A 2 B 2 . We need now to evaluate the three quantities A , 
B 2 and A - B, for which we will also use approximations based on the fact that 
L is soft while the other momenta are hard. We obtain first a rather simple, 
and by now very familiar [lqj, expression for A: 



A«Ao 2 * 2 



1 — cos ( 



2p 2 



2p 2 - 



2 2 
P p 



where 9 is the angle between the 3-vectors p and q, and with 



o 2 r 

Ml, + T-P(P + 9o) + U—p(p + q ) • 



(24) 



(25) 



7 I L is an analytic function of its arguments. The two possible choices for the square root 
of the complex number A lead to the same result. 
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Therefore, apart from the fact that the effective mass M^ s now gets an imaginary 
part coming from the width (and has been extended to hold for hard qo as well), 
this quantity is nothing but the denominator appearing in Eq. (39) of |l6[ ]. 

Since A ■ B prL 2 3> %/A and since the bremsstrahlung lives in the region 
where L 2 < 0, the difference of the two logarithms in Eq. ( p3| ) is just the 
discontinuity of the logarithm across its branch cut, which gives 



ln(A • B + VA) - ln(A • B - VA) w 2i7re(r) . 



(26) 



Therefore, in presence of a width T, the contribution of bremsstrahlung to the 
photon polarization tensor is 



Im%(Q) 



+ OC 



dp 



p 2 + (p + go) 2 



dip 
lo 



xDisc 



32tt 4 ql 

u 

l 3 dl [p T (l ,l)-p L (l ,l)] 
e(r) du 



[n F (p + q 

I, 



n P (p)] 



1 



(?)' 



1p* 



HI 



L 2 M 2 

— 5 1 

p 2 p 2 



1/2 



(27) 



where u = l — cos9. At this stage, this expression is formally similar to the one 
found for T = and reproduces known results in this limit. Indeed, when T = 0, 
then M 2 S is a real number, and taking the discontinuity just gives a factor of 
2 (do not forget the e(r) in the numerator). It is then possible to use the very 
same sequence of changes of variables as in p"6| ] to transform the integral over 
u, and write ImlP^Q) afl 



ImIP„(Q) 



2eV T 



-t-oo 



dp 



p 2 + (p + qo) 7 



m=T,L 



dw 



n F {p + qa) ~ n F (p)] 



(w + Rmf + {I m ) 2 



(29) 



where we denote: 



—L 2 



Re M: 



imir 



ReM e 2 ff 



- l ± 
I 



ImM c 2 ff 
ReMi 



Ren, 



Re Mi 



(30) 



8 We used the fact that the spectral functions p T L 
region: 



have a simple expression in the space-like 



XL) 



-2Im<7(L) 



{I? 



-Ren"7(L))2 



(28) 
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and where the function K comes from the integral over u and is defined as[] 



K(w,T) 



dy 



y + i/w 



VI -V (y + 4/w) 2 + {4T/ W y 



1 



4(a 2 
-2(3 



+ ft 2 
arctan 



{l + af+ft 2 



ft 



(l-a) 2 +ft 2 



arctan 



ft 



(32) 



with a + ift a square root^j of the complex number (4 + w)/w + AiT/w. This 
function K(w 7 T) is the generalization to the case of a non vanishing width of 
the factor y/w/(w + 4) tanh -1 y/w/(w + 4) appearing in Eq. (89) of In the 
limit of vanishing width (r — > 0), we recover the results of |^6[ . Eq. (|2^) cannot 
be further simplified analytically (except in some limiting cases), and must be 
evaluated numerically. 



5 Discussion 



5.1 Modification of the emission spectrum 

By inspecting our final expression given in Eq. (p9|), the first thing we notice 
is the very strong similarity with the same result in the absence of the width. 
Only the function K(w, F) contains the width, in the form of the dimensionless 
ratio r. When this ratio is small, the width has a negligeable effect while on 
the contrary if T 3> 1 then the width plays a dominant role. 

A first simple conclusion is obtained by noticing that the width arrives in 
M 2 S via the combination Tp(p + qq) / qq which becomes large at small go- There- 
fore, the effect of the width is more important at small qo. This property is 
illustrated by the plot of Fig. |^, which shows ImII M M (obtained numerically 
from Eq. (p9|)) as a function of the photon energy qo (the invariant mass Q 2 
is kept zero in this plot), for different values of the width T. On this plot, 
one can see that the q^ 1 behavior of bremsstrahlung at small go is modified in 
presence of a width. Instead of that, one reaches a plateauf^] when go - * 0. By 

9 One can go from the variable u to y by the following transformations 

2p 2 u = Af c 2 ff [w(l- z) + z- 1 - 1] , and y = 4(z - z 2 ) . (31) 



'Explicitly, we have: 



a,j3 



((4 + w)/w) 2 + (4r/wi) 2 ± ((4 + w)/w) 



(33) 



1 1 A similar behavior lias already been observed in a very different calculation by Weldon 
in |B0| (see Fig. 1 of [B0|). However, the resummations considered in this paper deal with 
the possibility to emit more than one photon, and affect the spectrum only below the scale 
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Figure 4: Effect of the width on the bremsstrahlung as a function of qo/T (for 
Q 2 = 0). The various curves correspond to different values of the width T. From 
top to bottom, the ratio TT/M^ (which corresponds to T/4 when qo ^> T) takes 
the values 10~ 6 , 1CT 4 , 1CT 2 and 1. 



comparing the second and third curves, we see that the value of qo at which 
we reach this plateau varies by a factor 10 2 when T varies by the same factor. 
This is a consequence of the fact that the width enters in the result via the ratio 
T/qo for small qo. Moreover, the value of the plateau is proportional to 1/r 
(this can be deduced from the plot since T increases by a factor 10 2 between 
two consecutive curves). This leads to the conjecture that the rate has a very 
simple r dependence in the region where the width is dominant. Additionally, 
when r increases, the extension of this plateau increases as well. In particular, 
for a large enough width, even the region of hard photons is modified. 

5.2 Region where the width is important 

If one wants to go beyond this qualitative statement, and get a better idea of 
the region (in the (q,qo) plane) where the width is important, one must study 
the ratio T. Indeed, the curve defined by the equation T = 1 is precisely the 
curve on which the effect of the width is comparable to the effect of and 
Q 2 (see Fig. ||) . Above this curve (region II of Fig. |), the width is an irrelevant 
parameter, and below this curve (region I - low invariant mass photons) we 
have the region where the width is the dominant collinear regulator. 

At first sight, the region where the width is important looks very similar to 
the region where higher loop corrections are IR-sensitive to the scale g 2 T and 
must be resummed, as found in Jl8| . We recall that in that work we considered 
the infrared structure of higher loop diagrams contributing to thermal photon 
production. It was shown that compensations between different cuts cancel 

Qo ~ a EM T, much smaller than the scale at which T ~ g 2 T\n(l/g) starts playing a role. The 
effect we are considering here appears at much larger photon energies (even including hard 
photons if r is large enough) , since it is related to resummations of gluons and thus involves 
the strong coupling constant ct s S> a EM . 
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q 

Figure 5: Boundary obtained from the condition T = 1. In region I, the width 
is the dominant regulator of collinear singularities. In addition, the shape of 
the emitted spectrum is modified by the width. In region II, the width is only 
a sub-dominant correction. 



all unscreened infrared divergences, and that the remaining terms are sensitive 
to gluon momenta down to a Q-dependent cutoff. In addition, we compared 
this cutoff with the scale n ~ g 2 T of the magnetic mass, and concluded that 
higher loop diagrams would be sensitive to the magnetic mass for a small enough 
photon invariant mass (see the figure 5 of ]l8fl). 

This similarity comes from the fact that line in the (q, go) plane where T 
starts to be important is given by the equation: 

2r ^ ReM^ (34) 

while the line on which the magnetic mass /i becomes important has the follow- 
ing equation: 

^ q j ReM \. (35) 

In other words, the two lines are defined by comparing a common momentum 
scale alternatively with the width of the quarks and with the magnetic mass. 
The physical interpretation of the momentum scale appearing in the right hand 
side of Eqs. ( |34| ) and ( |35|) will be given in the next section. The two boundaries 
are therefore very similar in QCD because T ~ g 2 T\n(l/g) and /i ~ g 2 T are 
not very different. But on the other hand, T and [i have very different mean- 
ings^ , and we may expect rather different interpretations for the two conditions 
Eqs. (34|) and (|3^). In the next section, we show that Eq. ( |34| ) is closely related 
to the LPM effect. 

There is also a more technical argument suggesting the difference of Eqs. ( [34) ) 
and d35|), which becomes clear by examining how the width T would appear in 

12 In particular, the fact that T and /i are very close is specific to QCD. In QED, one would 
have fi = and T ~ e 2 Tln(l/e). 
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perturbation theory. For that purpose, let us insert a self-energy correction 
on the leg of momentum R in the two-loop diagram of Fig. |^. Since we are 
dealing with the width perturbatively here, the quark propagators contain only 
the asymptotic mass M^. By a crude power counting, we can estimate that 
this insertion brings the following extra factor 13 



ir°T 



iTp(p + q Q ) 



R 2 - M 2 



ReM 2 s 



(36) 



We can already see that a self-energy insertion increases the strength of the 
potential collinear divergences by bringing an additional denominator R 2 — M^. 
On the other hand, this insertion does not modify the infrared properties of this 
diagram. By summing over the number of such insertions from to +oo, we 
get a factor 

(37) 



ReM e 2 ff 



ReM 2 s ®ir P (p + qo)/q () ' 

Therefore, the effect of such a resummation is to substitute Re M 2 S by Re M 2 S 
iTp(p + qo)/qo in the factor T 2 /M 2 S of collinear enhancement. This is precisely 
what has been observed in the more rigorous calculation of section [|. In ad- 
dition, this simple argument demonstrates clearly that the mode of action of 
width insertions is to affect the collinear sector, leaving unmodified the infrared 
sector. On the contrary, the possibility to have a sensitivity to p found in |jl8f 
is related to infrared singularities due exclusively to transverse gluonsp]. 



5.3 Limit of dominant width 

In region I, where the width becomes the dominant regulator, the ratio T is 
large, which enables us to make some additional approximations. In particular, 
we can perform very simply the integrals over w and y in Eq. (p9[). Indeed, 
we can first note that T sets the order of magnitude of the variable w. As a 
consequence, typical values of w are large, and we can neglect the corrections 
R T L and I T L in the denominator of Eq. (p9|), as well as 1/w in front of y. 
Therefore, this equation becomes 



i m n%(<2) 



2e 2 g 2 T 



dp 



p 2 + {p + qof 



m=T,L 



d.v 







+oo 



dw 
u>' 2 



[n F (p + q ) - n F (p)) 



dy 



o 



o 



V 1 -yy 2 + (4r/w) 5 



13 A self-energy correction of the type considered here modifies the real-part of the pole of the 
propagator by a mass-shift SM^ , negligible compared to arising from HTL corrections. 

14 In particular, longitudinal as well as transverse gluons contribute to T. There is no 
contradiction jsqth Jig ] which found only the transverse gluons to be important, since the 
statement of |18| was about infrared singularities. 
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2eV T f p 2 + (p + g ) 2 r , ( „ 3vrm 



ggy* 2 [n - (p + 90) -^ W]X 2R^iP (38) 



(i 



where m g ~ <?T is the gluon thermal mass coming from the prefactor of I T 
We see that this expression is of the form 



i{Q) ~ e y Y^VT) ' (39) 

where $ is a function independent of T giving the shape of the photon spectrum. 
Therefore, in the region where the width dominates, the spectrum scales as 
One notes that the combination THeM 2 s is proportional to Tp(p + qo)/qo, and 
therefore the shape of the spectrum is independent of Q 2 and M^. Furthermore, 
for qo <C T, it is independent of go (see Fig. ^). 

It may be interesting to compare in two extreme cases (qo soft and qo ^> T) 



the result of Eq. (39) with those obtained with a vanishing width. In the soft 



qo case, the result when r = was obtained in |lq : 

Imn%(Q)~eV — \, (40) 

where we have explicitly isolated the factor T 2 /M^ ~ l/<? 2 that comes from 
the collinear enhancement, while we have now: 

I mn %(Q)~eV-^- (41) 

The factor qo/F is in fact the new enhancement factor coming from T 2 /Im M 2 S . 
Therefore, the two results differ only by the nature of the factor of collinear 
enhancement. For hard photon bremsstrahlung, we had on the other hand [0 

ImnVQ)~eVr 2 ^, (42) 

which becomes now 

ImIPyQ)~eVT 2 ^ . (43) 

Again, the two expressions differ only by the collinear enhancement factor: 
T 2 /M^ - l/g 2 if T = instead of T 2 /ImM 2 s - T/T when the width is 
dominant. The numerical prefactors not written in Eqs. ( fl2| ) and (^) could 
however be quite different since they reflect a different physics. 

The r _1 scaling law in Eq. ([59j) has also a very simple physical interpretation 
related to the fact that is the mean free path of the quark in the medium. 
This result just tells us that the photon rate is proportional to the mean free 
path of the quark that emits the photon. In other words, quarks colliding very 
frequently do not have enough time to emit photons. This picture is at the basis 
of the Landau-Pomeranchuk-Migdal effect, that we discuss in the next section. 
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6 Connection with the LPM effect 



6.1 Generalities on the LPM effect 

The non-perturbative region found in the present paper can also be interpreted 
in a much more physical way, which seems to indicate that the region I in Fig. | 
is the region where the Landau-Pomeranchuk-Migdal ]3l], ||^, ||[ |34], |35| effect 
modifies photon production by a plasma. 

Let us first recall the condition for the LPM effect, by using a very heuristic 
argument valid for real photon production. In the bremsstrahlung process of 
photon production, the gluon exchanged during the scattering has a minimal 
momentum given by [ p3| 

go Ml 

2p(p + q ) ' ^ 
where p is the momentum of the quark. The inverse of this momentum transfer 
defines the "length" on which the photon is emitted, and for this reason is called 
the "coherence length" (denoted A co h in the following) for the production of this 
photon. It is also interpreted as the "formation time" of the photon. If A co h 
is much smaller that the typical distance between two consecutive scatterings 
(proportional to 1 /F) , then the photon rate is not affected by multiple scatter- 
ings. In other words, successive scatterings can be considered as independent. 
On the contrary, if the coherence length is larger than the mean free path, then 
successive scatterings are not independent anymore. This is the LPM effect. It 
is also possible to recast the previous condition as a comparison between the 
"formation time" and the mean free path of the quark in a very suggestive way: 
the photon must be produced before the quark scatters off another parton. 

6.2 LPM effect in thermal field theory 

Let us now show that the previous discussion arises automatically in the thermal 
field theory approach. In fact, there is in the calculation presented before a 
quantity very similar to the coherence length discussed in the semi-classical 
treatment of the LPM effect, namely, 

This value agrees with the one of H, |H, H when Q 2 = since it comes from 
kinematics, and generalizes it to the case where the invariant mass Q 2 is non 
vanishing. Additionally, this quantity was found in |Q to be the lower bound 
for the momentum of exchanged gluons. 

It is now possible to reformulate the condition T = 1 in more physical terms. 
Indeed, this condition can be rewritten as 

Amean ~ T ~ A C oh > (46) 

where A mea n is the mean free path of the quark in the plasma, and the region 
where the width dominates corresponds to the condition A mcan < A co h- At this 
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point, the connection with the LPM effect is obvious: the region we found to 
be non-perturbative due to the width of the quarks is also the region where the 
LPM effect matters. 

This discussion on the condition for the LPM effect to occur is in fact sum- 
marized in the expression for M 2 S introduced in Eq. ( p5| ) . This formula can be 
written very elegantly as: 



M\ = 2p(p + go) 



cff 



1 



(47) 



and the physical condition for the LPM effectfj£] to be relevant is mathemati- 
cally expressed by the dominance of ImAT^ over ReM e 2 ff . Conversely, when 
ReAf c 2 ff > ImM 2 ff , the perturbative approach is valid. As a side remark, let 
us note that the quantity M e 2 fi , which controls the physics of bremsstrahlung, 
combines in a very simple way three soft scales of the problem: Q 2 , and 
T. It is also possible to give a very simple physical interpretation of the ratio 
r which contains all the T dependence of the final result (Eq. (|29|)). Indeed, 
one can rewrite this ratio as T — A co h/A moan , which is nothing but the typical 
number of coherent scatterings necessary to produce a photon. 

The analogy with the standard treatment of the LPM effect is only partial 
though. Indeed, our estimate of the effect of the finite mean free path of the 
quarks is too crude (because of the constant width, and because of the missing 
vertex corrections) to be reliable from a quantitative perspective. The trend 
found here (suppression of the photon spectrum at small go) is however expected 
to be a solid prediction, and is in agreement with what is usually found from 
the LPM effect. 

The main difference with the usual calculation of the LPM effect comes from 
the treatment of scattering centers: one usually assumes a fast moving charged 
particle going through a medium of "cold" (static) scattering centers. In the case 
we are considering here, both the particle emitting the photon and the scattering 
centers are thermalized particles, of comparable momenta. The approximations 
made in J3l], 32, (3f| do not apply to this situation, and it is not clear 



whether a quantitative agreement is to be expected at all. In particular, the 
static approximation selects only Debye shielded longitudinal gluon exchanges, 
while the width (the resummation of which introduces multiple scatterings in 
our approach) of the thermalized quark receives contributions from transverse 



gluons as well. Contrary to |3l[ [32| [33|, |34j, |35| , transverse gluon exchanges are 
very important in thermal field theory, and matter for the LPM effect as well. 

There is a major difference between the LPM effect from multiple longitu- 
dinal gluon exchanges and multiple transverse gluon exchanges, which can be 

15 In [ferjl was investigated the influence of the finite mean free path of nucleons on axion 
(or neutrino pair) production by a supernova. The connection between the wid th V and the 
LPM effect is also mentioned in this paper. However, the condition given by Eq. (Ma) does not 
come out from the formalism of [B£j (it seems that the approximations made inp6[ for the 
bremsstrahlung do not enable one to track the coherence length, so that only T appears in the 
final result). What we have shown in the present section is that a careful calculation of the 
bremsstrahlung in thermal field theory generates automatically both terms of the comparison, 
through the effective mass given in Eq. (^). 
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readily seen when one compares the respective ranges of electric and magnetic 
fields to the mean free path of the quark. Indeed, since A mean » m^ 1 ~ (gT)^ 1 , 
the mean free path is much larger than the range of Debye screened electric 
fields. As a consequence, successive exchanges of longitudinal gluons are inde- 
pendent: they correspond to scatterings off different partons. On the contrary, 
we have /i -1 > A mcan , which implies that the range of magnetic fields can ex- 
tend beyond the mean free path of the quark. Therefore, successive exchanges 
of transverse gluons may correspond to scatterings off the same parton. This 
interpretation is supported by the infrared study performed in fisfl , where we 
found that only transverse gluons are causing trouble in the infrared sector, 
and that this problem occurs when A co h > ( see Eq. (j35|)). This condition 
means that the production of a single photon occurs on a distance larger than 
the correlation length of magnetic fields: this emission process is therefore able 
to probe the scale of the magnetic screening, and the rate is expected to become 
non perturbative. This qualitative difference of transverse gluon exchanges is 
also supported by the fact that the contribution of longitudinal gluons to T 
is perturbative (saturated at 1-loop), while the transverse contribution to T is 
non-perturbative |2^, |2^] . 



6.3 Comparison with the approach of Cleymans et al. 

The LPM effect has already been studied in the context of photon production by 
a quark-gluon plasma, although with a very different approach, in [j37| . In this 
paper, the study follows the semi-classical approach of p^| , to first deduce the 
photon rate from a single quark trajectory, and then average over the possible 
trajectories. Finally, the authors of |37i] manage to rewrite the rate as a function 
of a quantity F (9, r) which is related to the probability for the quark trajectory 
to undergo an angular deviation 9 after a time r. Additionally, they show that 
this object satisfies the following Fokker-Planck equation, which would be in 
our notations^: 

F = ^V 2 e F , (48) 



OF .q 



9> + M~ + % 
P % 



where is the mean square of scattering angle per unit length. 

In order to make the connection with our approach more transparent, we 
can factorize out the quark width in the following way (recalling the fact that 
r _1 gives the mean free path): 

(9 2 s ) = 16 9 2 T , (49) 

where 9 1 is the average square of the scattering angle per collision (up to a 
purely numerical factor) and where the prefactor 16 has been chosen for later 

16 In |J only soft photons are considered. As a consequence, in this comparison, we assume 
always qo >C p, and therefore do not distinguish p and p + qo. 
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convenience. We can now rewrite the equation satisfied by F in the following 
way 

dF 



r 6 2 



dr 



2p 



v 2 { 



ReM 



off 



F = 0. 



(50) 



We see that this equation is governed by a complex number whose real and 
imaginary part also appear in our approach. In particular, the condition for 
having LPM effect (discussed after Eq. (2.21) of ]37| ]) is the same as ours, given 

the fact that typical scattering angles satisfy 6 2 ~ ReM^g/p 2 . Note however 
that our calculation also includes the case of hard photon production, and shows 
that transverse gluon exchanges are equally important, both points which are 
not contained in Clcymans et al.'s approach. 



7 q*q annihilation 

7.1 Technical differences 

In section ^, when we performed the integral over the quark energy po using the 
S(P 2 — M~o) in Eq. |l9|, we considered only the positive energy solution, which 
corresponds to photon emission by the bremsstrahlung of a quark. The calcu- 
lation is very similar for the contribution of po — —uo p , with some peculiarities 
that we are going to highlight in this section. 

One of the most important differences is that, since we have qo > 0, this 
contribution contains in fact the sum of two processes (see Fig. ||) : the bremss- 
trahlung of an antiquark if uj p > qo, and the annihilation of an off-shell quark 
(a quark put off-shell by a scattering) with an antiquark if uj p < qo- The 
latter process has been studied in the case of zero width and real photons in 
[ p"7| and was indeed found to be a very important contribution to hard photon 
production. The purpose of this section is not to present a complete calculation 
of this process (which would require the discussion of a few other technical 
issues, especially in the region where Q 2 > 4M^), but only to reproduce in this 
case the discussion of section ||, in order to have a qualitative picture of the 
influence of the width on this process. 

The calculation of the po < contribution can be performed by the same 
method as the po > one. In particular, the angular integral over c?f2; generates 
an effective mass M^g, the expression of which is 

M' c 2 = Ml + %p{p - q ) + 4i-p(p - q ) • (51) 

In addition to this change, the angular integral over d£l p is now dominated by 
values of 9 around it (because 1— cos 9 in the expression of A becomes 1+cos 6>)p|. 
A practical consequence of this is that r w \p — qo\ in the collinear limit. All 

17 This difference is also obvious when one compares the process on the left of Fig. ^ with 
the two processes on the right of the same figure, which we are dealing with now. 
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the subsequent steps can be reproduced and one is lead to the following final 
formula: 

+00 

ImII% (Q) fa 2 / dp [n F (q -p)~n F (-p)j 

" <?o J * 


V- fdx T, I™ K(W,T) 

x V / — dw » = , 52 

where the notations are the same as in Eq. (^), with M' e $ in place of M^ s . 

Another very important technical difference appears when one performs the 
integral over the quark momentum p. Indeed, the range of this integral is 
controlled by the statistical factors n F (r ) — n F (p ) s» n F (p) — n F {p — q ). This 
factor is of order 1 in a domain going from p = to p ~ Max(q ,T). In other 
words, if qo is soft, this range is the usual [0, T] interval. On the contrary, if 
q ^> T, the integral extends to p ~ q ^> T. A practical consequence is that for 
qo <C T, the contribution of po < comes mainly from the bremsstrahlung of 
an antiquark (since most of dp integral satisfies qo < p ~ T), while for q Q ^> T 
it is dominated by the q*q annihilation (since we then have p < qo). In the 
intermediate region qo ~ T, this contribution is a mixture of both processes. 



7.2 Modifications due to the width 

Again, we see that the term depending on T in M'Jg comes with T/qo, indicating 
that the width will modify significantly the rate when qg is soft. But, in addition, 
when qo ^> T, p can itself be of order qo, so that the imaginary part of Af^J is of 
order Tqo, which can be much larger than at very large qo- It means that 
the width also changes dramatically the spectrum of very hard real photons (in 
the case of bremsstrahlung, the width modifies only marginally the rate of hard 
photons, since the imaginary part of M^ s is of the same order of magnitude as 
if qo ^> T and T ~ g 2 T) coming from the q*q annihilation. 



This has been checked by evaluating numerically Eq. (52), and the results are 
displayed in Fig. ||, where we plot the imaginary part of the photon polarization 
tensor as a function of qo/T (with the same values of temperature and coupling 
constant as in Fig. |i|) for a set of values of T identical to the one chosen in 
Fig. H From this plot, it is obvious that the regions qo <§C T and qo 3> T are 
very different. In fact, in the region qo ^ T, we obtain results that exactly match 
those of Fig. |J indeed, in this region, the po < contribution is dominated by 
the bremsstrahlung of an antiquark, which as expected contributes equally as 
the bremsstrahlung of a quark. 

The new feature coming with po < appears in the hard photon region 
(qo 3> T), where now the q*q annihilation contributes. At very small T (upper 
curve), we recover the result of fTij ] according to which the process q*q —> 7 
gives a contribution to ImlF'^Q) that increases like qo- Then, as F increases, 
we see a saturation at large enough qo, to a value that behaves like This 
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Figure 6: Effect of the width in the region po < as a function of qo/T (for 
' > J — 0). Each curve corresponds to a different value of the width T. The set 



of values taken for T is the same as in Fig. ||. 



confirms the above qualitative statement based on M'J^. In this case also, it is 
possible to obtain a formula like Eq. ([39]) showing explicitly the V scaling law 
when r is the dominant collinear regulator. 

The contribution of the process q*q — *■ 7, first evaluated at leading order in 
[ p"7[ , has been found to be phenomenologically important for the production of 
direct photons in the 2 — 10 GeV range, because it dominates previous estimates 
by a factor of order 5 p8[ . The inclusion of this contribution in heavy nuclei 
collisions simulations ]39[ has lead to a good agreement with the measured rates 
from the WA98 experiment Q. These rates did not include the LPM 
suppression advocated here, but one must realize that the energy range where 
thermal photons are relevant is also the region of minimum sensitivity to the 
width r (around the minimum of the curves in figure |]) , since for a temperature 
of a few hundred MeV, the GeV range corresponds to qo/T ~ 10. It is too early 
to be more quantitative here given the fact that vertex corrections 0] have been 
completely disregarded in the present work, but having a precise prediction for 
this process would definitely be of important phenomenological interest. 



7.3 LPM effect in q*q annihilation 

When po < 0, the condition under which the width is the dominant collinear 
regulator can be rewritten as 



2r > 



go BcMg 
2 p\qo ~p\ 



(53) 



and the quantity in the right hand side is the minimal value Z m i n for the mo- 
mentum L of the gluon exchanged in the scattering. Therefore, its inverse is 
also the coherence length for the emission process, and the above condition is 
nothing but the criterion for having LPM effect. 



22 



As a consequence, we are lead to the conclusion that the LPM effect plays a 
role not only for the photon production by bremsstrahlung at low energy, but 
also for the q*q annihilation at high energy. This conclusion was also recently 
obtained in [[l3| by a semi-classical method following the original method of 
Migdal [[33|. The authors of Q in fact studied the inverse of the process 
q*q — > 7, namely the production of a fermion pair out of a photon (this process 
is made possible by the fact that one of the fermions is produced off-shell and 
then scatters in the medium). They found that the LPM effect is important 
for the fragmentation of very high energy photons. This is equivalent to our 
observation that the LPM effect modifies q*q — > 7 for large go- 

8 Conclusions 

In this paper, we have studied the effect of the finite mean free path of quarks 
on photon production by a quark gluon plasma. The main result is that such a 
correction may have important effects by affecting the collinear sector. 

We have first considered the effect of a quark width in the simplest 1-loop 
diagram, and shown that potentially large 1-loop contributions cancel when a 
more complete summation is performed. This cancellation is not new in fact, 
and has been noticed elsewhere. It has connections with the fact that there are 
no hard thermal loops for vertices with 2 photons and any number of gluons. 

Two-loop contributions beyond the eikonal approximation escape this can- 
cellation, and lead to a result which exhibits features of the LPM effect. It is 
indeed possible to interpret the region where the width of the quark is impor- 
tant as the region where the formation length of the photon is larger than the 
mean free path of the quark. Our study also shows that the LPM effect modifies 
the emission spectrum of very hard photons. Despite its nice interpretation, our 
present result is not complete because it does not consider the vertex corrections 
that should accompany the resummation of a width. 

Multiple scatterings due to longitudinal gluon exchanges are independent, 
as in the semi-classical treatment of the LPM effect. On the contrary, the trans- 
verse sector displays features that are qualitatively different from what is usually 
accepted in the semi-classical treatment (which just disregards transverse gluon 
exchanges). Indeed, we find that multiple transverse gluon exchanges are very 
important also, and that they are not independent due to the long range of 
magnetic fields. 

The various scales in the problem are summarized in Fig. ^. Three of these 
scales are intrinsic to the quark gluon plasma: the Debye screening length m~ , 
the mean free path of the quarks A mcan ~ and the magnetic screening 
length The fourth scale, the coherence length A co h (or formation length 

of the photon), depends on the energy and invariant mass of the photon one 
wants to observe, and should be compared to the first three. This leads to three 
zones which have simple physical interpretations. In addition, the nature (and 
complexity) of the dominant higher loop corrections depend on the coherence 
length. 
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Figure 7: Summary of the various scales in the problem, and comparison with 
the coherence length. In region A: < A co h < A moan . Region B: A mca n < A co h < 
\i~ x . Region C: < X co h- Each picture shows the new class of diagrams that 
one must consider when going to larger and larger coherence lengths (in addition 
to the diagrams already considered for smaller coherence lengths). A boldface 
quark line denotes the resummation of a width in addition to the asymptotic 
thermal mass. 



• Region A: < A co h < A moa n, and photon production is dominated by single 
scatterings. The large scale structure of magnetic fields is irrelevant. The only 
relevant diagram is the 2-loop diagram (with HTL resummed quarks) already 
considered in |l6| . The contributions of transverse and longitudinal gluons are 
comparable. 

• Region B: A mcan < A co h < H , which implies that photon production is af- 
fected by multiple scatterings (LPM effect). In this region, the emission process 
is not yet sensitive to the magnetic screening. Indeed, according to this 
sensitivity comes in when A co h > M • It is sufficient to resum a width (satu- 
rated at Tlooppf) on the quarks, and to consider vertex corrections involving 
both longitudinal and transverse gluons. So far, we have said nothing about the 
vertex corrections that should come with self-energy insertions. Nevertheless, it 
is reasonable to assume from Ward identities that 1-loop self-energy insertions 
talk to ladder corrections^ We can note also that the discussion of this region 
is somewhat academic, because its extension is only proportional to ln(l/g). 
Additionally, the new topologies of the region C are at most suppressed by 
powers of ln(l/<7) _1 if evaluated in region B. 

• Region C: /i -1 < A co h. The LPM effect still modifies photon production. In 
addition, the emission of a photon lasts long enough for the process to be sensi- 
tive to the magnetic screening. This is the physical meaning of the result of Jl8j . 
Diagrams with an arbitrary number of ultrasoft transverse gluons, connected in 

18 When the IR cutoff is ^> /x, corrections to Y coming from topologies beyond 1-loop 

are suppressed by g 2 T/\^ h < 1. — — 

19 The mechanism which makes these ladder corrections important is the same as in |24[ |25| , 
and has nothing to do with infrared singularities. The arguments of Jis) cannot exclude 
them, even if the infrared cutoff is X~ o]l S> fi. This is why these corrections are due to both 
longitudinal and transverse gluons. 
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all the possible ways, must be resummed. In addition, the width included on 
the quarks contains non perturbative contributions due to transverse gluons 
[ p2| p3[ , and the effective gluon vertices must be corrected to hold for ultrasoft 
momenta p4j . 

A full study (including vertex corrections) of the corrections to photon pro- 
duction due to longitudinal gluons seems within the reach of perturbation the- 
ory, but the situation is very different in the transverse sector which seems far 
beyond the possibilities of perturbative methods. In this respect, photon pro- 
duction is not very different from the calculation of the quark damping rate 
[|H New tools, like functional methods, transport equations, and eventu- 
ally lattice techiques, are presumably the way to go in this area. In particular, 
the picture emerging from Fig. |?| suggests an analogy with the successive resum- 
mations of in which the scales T, gT and g 2 Tln(l/g) are successively 

integrated out. Indeed, going to larger and larger coherence lengths requires the 
inclusion of more and more complicated topologies, a procedure which amounts 
to integrate out degrees of freedom encountered at smaller length scales. 
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A Separation of the various cuts 

Equation (^) includes automatically the sum of the cuts (c) and (d) of Fig. |[ 
To see it explicitly, let us first denote 



F(r) 
G(r) 



At 2 
2p2 



e(r) 



JZ. 

2p 2 



p2 p2 



1/2 



(54) 



Then, the discontinuity on the last line of Eq. ( p7| ) equals: 

DiscF(r)G(r) = F(T)G(T) - F(-T)G(-T) 

[G(r)+G(-r)l , . \F(T) + F(-T)] . , , . 

= o DiscF(r) + 1 V ; V ^ DiscG(r) . (55) 



In the last line, the first term corresponds to cut (d) and the second term gives 
cut (c). The fact that both cuts contribute to the photon rate when f ^ was 
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to be expected. Indeed, the cut (d) which corresponds to the process q — ► 97 
is kinematically forbidden if T = 0, but is now allowed by the very fact that 
the width T reflects the collisions of the quasi-quark (and the channel q — * q-f 
is allowed in the presence of a medium). 

The integrals to be calculated are more complicated if one wants the two 
cuts separately. But it turns out that the difference between the two cuts is also 
very simple. This difference is given by: 

(c) - ( d) oc \m±iti)i DiscG(r) _ m+ga DiscF(r) 

= F(-T)G(T) - F(T)G{-T) . (56) 

Therefore, in order to get (c) — (d), one has to start from Eq. ( f27j) in which the 
last line is substituted by Disc F(— F)G(r). The same transformations can be 
applied to this integral, and one finally obtains an equation similar to Eq. (^9|), 
but in which the function K(w,T) is replaced by 

L(t0) f) s l / y + A/w x . (57) 

This new integral is also elementary and can be obtained in closed form if 
needed. 

It is obvious that in the limit of zero width (T — > 0), the two integrals 
K(w,T) and L(w 7 T) become equal. The only way for this to happen is that 
the contribution of cut (d) goes to zero. Therefore, we recover the fact that the 
process q — » qj disappears when T = 0. 
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